Solid state physics of impact crater formation: further considerations 
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Impact craters exist on solid surface planets, their satellites and 
many asteroids. The aim of this paper is to propose a theoretical ex- 
pression for the product pr^u^, where the three symbols denote the 
mass density,radius and speed of the impactor. The expression is de- 
rived using well known results of solid state physics, and it can be used 
in estimating parameters of impactors which have led to formation of 
craters on various solid bodies in the Solar System]^ 
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1 Introduction 

Craters of various sizes have been observed on the terrestrial planets,their 
satelhtes and the major part of the asteroids. The study of these craters, and 
the resulting constraints on the related impactors has become a separate 
field of research in planetary science, see for example [8]. At the time of this 
writing, the latest example testifying about the importance of impacts into 
the Earth and their consequences is the small asteroid which entered the 
atmosphere over the city of Chelyabinsk in Russia on February 15, 2013. 
One of the fundamental questions concerning the impactors is what can be 
concluded about them by combining astronomical data with results of solid 
state physics. 

Recent theoretical work for example [1] has shown that the application 
of elementary principles of solid state physics to this problem gives physi- 
cally plausible results in reasonable agreement with those obtained by celes- 
tial mechanics. It was assumed there that the material of the target was a 
crystal lattice, and the calculations were performed per unit volume. It was 
assumed in that paper that the condition for the formation of a crater is 
that the kinetic energy of a unit volume of the impactor has to be equal 
to the internal energy of a unit volume of the material of the target. The 
result was an expression for the minimal speed which an impactor of given 
parameters must have when hitting a target with a predefined set of pa- 
rameters, in order to form a crater. Obviously,a certain part of surfaces of 
planets, satellites, asteroids,., in which impacts occur are granular. For recent 
theoretical work on "granular impact" see, for example, [3] or [11]. 

The aim of the calculation reported here is to take into account the di- 
mensional effects - to consider both the impactor and the crater it forms as 
objects of finite dimensions. The novelty of the approach discussed in the 
next section, compared to existing work such as [4], is the generality - it is 
based on principles of solid state physics, and it can be applied to any solid 
material. Formation of impact craters is here discussed from the viewpoint 
of pure solid state physics as the following analogous problem: what kinetic 
energy of an impactor is needed to produce a hole of given dimensions in a 
material with a predefined set of parameters? The calculations will be per- 
formed under the assumptions that the material of the target is a crystal,that 
one of the five usual bonding types exists in it and that in the impact the 
target does not heat up to its melting temperature, so that solid state physics 
can be applied. For recent work on the problem of heating in impacts see [2] . 



The expression for the cohesion energy used in the present work is one of the 
theoretically possible expressions. More details are available, for example, in 

[7]. 

2 Calculations 

The physical keypoint of the calculation to be discussed in this section is the 
idea that the kinetic energy of the impactor must be greater than or equal 
to the internal energy of some volume, V2 oi the material of the target. The 
kinetic energy of the impactor of mass rrii and speed Vi is obviously: 

Ek = -mivf (1) 

The internal energy of a volume V2 of the target material is equal to the sum 
of the following three "contributions" : the cohesion energy ii^^the thermal 
energy Et and the energy Eh{T) required for heating the specimen by an 
amount AT in the moment of impact: 

Ei = Ec + Et + Eh{T) (2) 

The condition for the formation of a crater as a consequence of impact of a 
projectile into a target is, in general terms, given by 

Ek = Ei (3) 

In order to ensure the stability of matter, the interatomic interaction energy 
must have the form 

E = ^^^ (4) 

where R is the inetartomic distance, the first term on the right hand side is 
attractive and the second one repulsive.lt can be shown [5] that the cohesion 
energy is given by 

Ec = -^^^ (5) 

pn 

where Bq = —V{dP/dV)T = p{dP/dp)T is the bulk modulus of a material 
and Vim the volume per pair [5]. The speed of sound waves in a material with 
pressure P and density p is given by 
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The analytical form of the function Et depends on the temperature of the 
solid. It can be shown (for example [7]) that for temperatures 

, ^ ^w , , 
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{h is Planck's constant divided by 27r, u is the mean speed of sound in the solid 
and a is the lattice constant), one should use the "high temperature" limit 
for the thermal energy. Inserting for example u = 5km/ s and a = 5 x 10~^m 
it follows that Tq » 77 K. 

The thermal energy Et is given by 

Et = 3NukBTD3{x) (9) 

where x = To/T, T^, is the Debye temperature,// is the number of particles 
in the elementary crystal cell,iV is the number of elementary crystal cells in 
the specimen, and ks is Boltzmann's constant. The symbol D^{x) denotes 
the n = 3 case of the Debye function, which is given by 

^ ^ /- tVt _ ,„.„.!_ X ^ B2kX^' 

"^^^ x« io Exp[t] - 1 ^""/^ ^^ ^n 2(n + 1) ^ f^, {2k + n){2k)V 

(10) 
and i?2fc denote Bernoulli's numbers. Integrating and retaining the terms up 
to and including second order leads to: 

Et = SNuksTil - ^^ + (^)(^)1 (H) 

The specific heat Cy is given by 

Cv = ^ = 3NukB{Ds{TD/T) - (Td/T)D',{Td/T)) (12) 

which leads to the following approximate result: 

Cv = 'iNukB\l-(—)(— ? + (—)(— f] (13) 
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The energy Eh{T) is obviously given by Eh = Cy(Ti — T), where Ti is 
the temperature to which the material of the target heats up in the moment 
of impact and T is the initial temperature. 

One of the problems consists in knowing the volume V2, which is, in first 
approximation,equal to the volume of the crater. It is known from theoretical 
work [4] and references given there that the form of a crater is determined by 
a combination of "gravity scaling" and "strength scaling", where the term 
"strength" reffers to the material strength of the target. Explicite expressions 
for the form of craters exist in [4] but they are applicable to 4 material types. 
In order to simplify somewhat the calculations, and get slightly more general 
results, it was assumed in the present work that craters have the shape of 
a half of a rotating elipsoid,with distinct semi axes, denoted by a,b and c. 
Physically,a and b denote the semi axes of the "opening" of the crater, while 
c is the depth. In that approximation,the volume of the crater is obviously 
given by 

V2 = -Txabc (14) 

o 

Inserting the definition of the mass density (pi = y^) in eq.(l) and assuming 
further that the impactor has the form of a sphere of radius ri,the following 
expression for the kinetic energy of the impactor is obtained: 

Ek = -mivf = - X -Trrfpivf = -Trpirfvf (15) 

Inserting eqs.(7), (11), (13) and (15) into (3) one gets the final form of the 
energy condition which has to be satisfied for the formation of an impact 
crater: 

^560^^ T^ 420T3Ti npukeTi^ 3 ^^ ^ ' 

The number A^ is equal to the ratio of the volume of the crater, approximated 
by eq.(14), and the volume of the elementary crystal cell,fe: A^ = V/Vg. 
Equation (16) groups known or measurable parameters on the left hand side, 
while on the right hand side it contains parameters of the impactor. Applying 
this equation demands the knowledge of the equation of state (EOS) of the 
material of the target. Generally speaking,regardless of its detailed form, an 



EOS can be expressed in the following analytical form 

p{p) = T.^^{-r (17) 

where a^ are some coefficients, po is the density at some pressure Pq and p is 
the density at pressure P. It follows that 

(|^)=f:(^+i)«m(fr (18) 

A well known example of the EOS of a solid is the Birch-Murnaghan EOS 
|12] 
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{l + (3/4)(5^-4)[(p/po)^/^-l]} (19) 

The symbols po and p denote the mass density of a specimen under consid- 
eration at the initial value of the pressure Pq and at some arbitrary value 
P. 

More refined results could be obtained by using analytic approximations 
to the Helmholtz free energy, from which all the thermodynamic potentials 
can be derived. This is, for example, the approach used in the ANEOS EOS 
[13]. However, ANEOS uses a file of up to 40 parameters for the characteri- 
zation of materials, which complicates its applications. 

3 A test example 

As a test of the applicability of the procedure discussed in the preceeding 
paragraph, it was applied to the Barringer crater in Arizona. One of the min- 
erals found in and around the crater is Forsterite - Mg2Si04. The dimensions 
of the crater are a = b = 1.186A;?t2 and c = O.nkm. Take that the expo- 
nents in the expression for the interparticle potential are p = 1 and n = 9. 
Assume furher that V = 5A;m/s,which is the mean value of the measured 
velocity of seismic waves in the Earth, that the temperature before impact is 
T = 300K and that the target heats up to Ti = 550K at the point of impact. 
This is safely far from the melting temperature of Forsterite (which is around 



Tm = 1900C) so that solid state physics can be apphed. All the parameters 



of Forsterite exist in the literature, such as [10] or |http://webniineral.com 



Taking the necessary values, and assuming that the impactor had a radius of 
ri = 65m and density of pi = 8500kg/m^ one finally gets that vi = Alkm/s. 
This value is larger than existing results, such as for example [9] or [6]. The 
discrepancy can be traced to the assumption, made in the calculation re- 
ported here, that the target material is only Forsterite. In reality, this is 
certainly not the case. Taking that 10 percent of the target is Forsterite, and 
that the object had a radius of ri = 60?7i would give vi = 15km/ s , which is 
much closer to the results of celestial mechanics. 

4 Discussion 

In this paper we have presented a simple procedure which gives the possi- 
bility of estimating the value of the product of the density and radius of the 
impactor and the speed of impact: pirfv"^ in terms of various material pa- 
rameters of the target and the impactor. The expression for this product has 
been derived using basic principles of solid state physics, without any spe- 
cial assumption(s) about the materials. In (at least some of the) terrestrial 
applications the density of the impactor and the inclination of the trajectory 
can (in principle) be measured or estimated, which gives the possibility of 
estimating the impact speed if the value of ri can somehow be estimated. 

The shape of the craters has been approximated as a half of a rotating 
elipsoid,regardless of the composition of the material of the target. As a 
first approximation, the speed of the seismic waves has been "put in" in- 
stead of being calculated from the equation of state of the material of the 
target. Calculating this speed would have demanded the precise knowledge 
of the chemical composition of the target materials, and of Bq and Bq for 
them. The impact of a projectile into a target leads to heating,and possi- 
bly melting and even vaporisation of the target around the impact point. 
The thermodynamic result of an impact depends on the heat capacity of 
the material of the target and on the kinetic energy of the impactor. The 
heat capacity can be measured, or theoretically estimated, assuming prior 
knowledge of the chemical composition of the target. Knowledge of the heat 
capacity is vital for estimates of the temperature to which the target heats in 
the impact. Therefore, the question is which impacts lead to melting and va- 
porisation of the material of the target, and which just provoke heating of the 



target. Even if the material of the target gets partially melted and as such 
flows away, the dimensions of the crater formed by the impact are determined 
by the volume of the material pushed away in the impact. This volume is, in 
turn, determined by the ratio of the kinetic energy of the impactor to the 
internal energy of some volume of the target. Some details of this problem 
have been discussed in [2]. 

The general conclusion of work reported here is that by using simple 
well known results of solid state physics and a given appoximation of the 
shape of impact craters, it becomes possible to estimate the value of the 
product pirf{vcos9Y. This in turn can be used to draw conclusions about 
the impactors which made various craters on the solid surfaces of objects in 
the planetary system. 
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